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Abstract 

The properties of diluted symmetric A-B diblock copolymers at the interface be- 
tween A and B homopolymer phases are studied by means of Monte Carlo (MC) simu- 
lations of the bond fluctuation model. We calculate segment density profiles as well as 
orientational properties of segments, of A and B blocks, and of the whole chain. Our 
data support the picture of oriented "dumbbells", which consist of mildly perturbed 
A and B Gaussian coils. The results are compared to a self consistent field theory 
(SCFT) for single copolymer chains at a homopolymer interface. We also discuss the 
number of interaction contacts between monomers, which provide a measure for the 
"active surface" of copolymers or homopolymers close to the interface. 
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1 Introduction 



Blending of polymeric substances is a straightforward and inexpensive way of 
creating new materials with improved mechanical properties^]. However, long 
polymers of different type A and B are often immiscible already at high tem- 
peratures, since the total energy of the (usually repulsive) relative interaction 
is proportional to the total number of monomers and cannot be balanced by 
the entropy of mixing, which is proportional to the number of polymers in 
the mixture || In order to overcome this problem, block copolymers con- 
taining both types of monomers can be used as effective compatibilizers^, ||. 
Being partly compatible with both the A and B rich phase, they tend to ag- 
gregate at interfaces, where they reduce the number of direct contacts between 
A and B homopolymers, thereby reducing the interfacial tension 0. Conse- 
quently the total area of interfaces increases, and the immiscible components 
may get finely dispersed in the mixture. Furthermore, copolymers improve 
the mechanical properties of such interfaces: Due to entanglement between 
homopolymers and copolymers, they increase the adhesive attraction and the 
fracture toughness [0, || |]. At high enough copolymer concentrations, addi- 
tional copolymer rich phases emerge which may display a diversity of structures 
ordered on a mesoscopic scale [1C], [TT| . 

The simplest possible copolymers are diblocks, which consist of a block of 
A monomers connected to a block of B monomers. From a thermodynamic 
point of view, their effect on interfaces can be described as follows: Copolymers 
act as amphiphiles in the homopolymer mixture [|T2[. Increasing the copolymer 
concentration causes the interfacial tension to decrease monotonically, until 
this process is terminated by the formation of a third, copolymer rich phase, 
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e.g. an ordered lamellar phase or microemulsion. The amphiphilic strength 
of a copolymer, i.e., the maximum reduction of the interfacial tension it can 
achieve, increases with the length of the A and B blocks relative to the size 
of the homopolymers. Mean field theories predict that the interfacial tension 
can be driven to zero for long copolymers, which implies that the copolymer 
rich lamellar phase evolves into the two-phase region via a continuous un- 
binding transition]!^, H]. Fluctuations push the transition to first order in 
real systems fT5[|. One can argue that the low interfacial tensions at the pres- 
ence of copolymers result from the affinity of the system to such an unbinding 
transition [|T6|j. 

This macroscopic discussion however does not shed light on the microscopic 
mechanisms, why and how copolymers work as compatibilizer or amphiphiles. 
In a simple microscopic picture, the A block of copolymers at interfaces prefer- 
ably stick into the A-rich phase, the B block into the U-rich phase, and the 
molecules as a whole act as reinforcing rods J?]. The conformations of the 
copolymers determine the properties of the interface. A more detailed micro- 
scopic description of copolymer properties at interfaces has been developed 



by Leibler|T7]] and refined by Semenov [p!8| . The analysis is based on the as- 
sumption, that the junction points between A and B blocks are confined to a 
narrow region, which is much smaller than the width of the copolymer layer 
at the interface. The presence of copolymers at the interface gives rise to two 
free energy contributions: the entropy of mixing, divided into a translational 
term and a swelling term, and the elastic energy of stretching of the copolymer 
blocks. A scenario emerges which distinguishes between four different regimes: 
In the dilute regime, copolymers aggregate at the interface, but do not yet 
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overlap. The A and B blocks are described as weakly perturbed coils, and 
the free energy is dominated by the mixing energy. When the mean interchain 
distance gets of the order of the block radii of gyration, copolymers start to 
stretch and form a "wet brush" . At even higher copolymer concentrations, the 
"dry brush" regime is entered, where homopolymers do not penetrate into the 
interfacial region; A and B blocks have the conformation of stretched coils and 
the free energy is dominated by the elastic energy of stretching. Finally, in 
the saturated regime, the density of copolymers at the interface is close to one 
and the interfacial segregation of copolymers competes with the formation of 



micelles in the bulk 18 



Experimental studies of copolymers at homopolymer interfaces have been 
carried out by numerous groups, mostly using neutron reflectivity or forward 
recoil spectroscopy |19|-p5|. By deuterating individual parts of the homopoly- 
mers or the copolymers, the excess of copolymers at the interface can be mea- 
sured as well as distributions of homopolymer segments, copolymer segments 
or even selected copolymer segments, like the junction point between A and B 
or the end segment [|25| , |26| , [27| . These studies have provided detailed insight 
into the microscopic structure of such interfaces. 

Simple Leibler type theories are already quite successful in reproducing 
many of the experimental results flT8|, In order to reach quantitative agree- 
ment, however, the Flory Huggins parameter x has to be treated as adjustable, 
molecular weight dependent parameter. Less transparent, but more accurate 
mean field approaches are the self consistent field theories |28|, p9|j or the den- 
sity functional theories [50, |3l| . The self consistent field theory has first been 



applied to copolymer- homopolymer interfaces by Noolandi and Hong[Q and 
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adjusted to the case of a blend without solvent by Shull and Kramer p3|, |3~3 



It has been shown to be quantitatively successful in predicting the correct 
copolymer excess at the interface, using a x-parameter which is taken from 



independent bulk measurements |25| . The calculated width of the segment 
interface is somewhat too low, but the discrepancies can be understood quan- 
titatively if broadening due to capillary waves is accounted for[^2[ |25fl . Such a 
remarkable success of a mean field theory is characteristic for polymeric sub- 
stances, in which high molecular weight polymers interact with a high number 
of other polymers^, 35]. In addition to reproducing experimental data, the 
self consistent field theory also has the advantage of yielding further structural 
information, e.g. on chain conformations, monomer orientations etc.[]36], |15[] , 
which may be hard to access experimentally. Unfortunately, the theory also 
has some serious drawbacks. In particular, the usual treatment of polymers as 
Gaussian random walks is questionable for chains of the minority component, 
e.g. A in the B rich phase, and generally on length scales smaller than the 
screening length of the excluded volume [|45| . 



Computer simulations provide another way of obtaining additional informa- 
tion on the microscopic structure of interfaces. Simulations of inhomogeneous 
polymeric systems are computationally extremely demanding, and therefore 
rare. Minchau et al and Fried and Binder |38|, |39| have investigated the 
phase behavior of pure copolymer systems. Pan et al studied microphase struc- 
tures in systems of short copolymers, which are swollen by a small amount 
(volume fraction 10 %) of longer homopolymers[ Wang and Mattice have 
studied the adsorption of self avoiding copolymers at a stationary, sharp in- 
terface, modelled by an external field with a sharp kink[(|T|]. A similar study 
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has been performed for a random copolymer by Peng et al | i2| . More detailed 
simulations have been presented by Cifra[fi3|]. 

In this work, we present a study of copolymers at the interface between 
homopolymer phases, where both homopolymers and copolymers are treated 
in microscopic detail. Pure homopolymer interfaces in immiscible A/B blends 
have been analyzed previously [(44]] and compared to the predictions of a self 
consistent field theory |45|. Here, we consider the effect of adding a small 
number of symmetric diblock copolymers to such a "known" homopolymer in- 
terface. We restrict ourselves to the diluted case, where the copolymer coils 
do almost not overlap, and where the static structure of the interface is essen- 
tially that of a pure homopolymer interface. The results are compared to self 
consistent field theory calculations for a single copolymer in a homopolymer 
interface. Our paper is organized as follows: The next section describes the 
simulation model and method, and ends with a few comments on the self con- 
sistent field theory. The results are presented in section three. In particular, 
we discuss segment density profiles, chain and bond orientations, and profiles 
for the number of interacting contacts between monomers. We summarize and 
conclude in the last section. 

2 The Bond Fluctuation Model 

With the presently available computational resources, molecular modelling of 
the phase behavior in polymer blends in atomistic detail is far beyond feasibil- 
ity. Fortunately, many important features of such systems are already apparent 



in coarse grained modelsp6[, such as the bond fluctuation model on a cubic 



lattice |I7| . The latter models polymers as chains of N effective monomers, 
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which occupy each a cube of 8 neighboring sites and are connected by bond 
vectors of length 2, a/5, Vd, 3, or vTO in units of the lattice spacing a . One 
such cube represents a group of n « 3 — 5 chemical monomers. Hence a total 
chain length of 32, as has been used here, corresponds to a degree of polymer- 
ization of approx. 100-160 in a real polymer. At volume fraction <p = 0.5 Oq 3 or 
monomer density p = 1/16 Oq 3 , the model reproduces many important proper- 
ties of dense polymer melts, e.g. single chain configurations show almost ideal 
Gaussian chain statistics, the single chain structure factor follows a Debye 
function, and the collective scattering function has the experimental form . 

The relative repulsion between monomers of type A and B is modelled 
by introducing symmetric energy parameters Eaa = £bb = —£ab = —ksTe, 
which describe the pairwise interaction between monomers at distances of less 
than \/6«o- At e = 0.1, the interactions between monomers are dominated by 
the effect of excluded volume, and the mixture can be described as a weakly 
perturbed athermal melt; in particular, the equation of state and the compress- 



ibility are almost not affected by the presence of the interactions [|49|. From ex- 
tensive previous study of this model, the relation of these model parameters to 
commonly used parameters in polymer theories is well known. At chain length 



iV = 32, the statistical segment length b is given by b = 3.05 aofH]; the radius 
of gyration is R g = yN/Q b ~ 7 a , the compressibility is k B Tn = 3.9 aol49|j, 



and the Flory Huggins parameter \ can be calculated using x — 2££, where 
z = 2.65 is the effective coordination number in the bulk, i.e., the average 
number of interchain contacts of a monomer [|5T| . 

The interfacial properties were studied in a L x D x L geometry at system 
size D = 64 and L = 512. The dimensions of the system were chosen such 
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that the width D of the slab is much larger than the gyration radius R g , i.e., 

almost ten times as large. The boundary conditions are periodic in x and z 

direction and "antiperiodic" in the y direction, i.e., A-chain parts leaving the 

right part of the simulation box reenter it on the left side as S-chain parts 

and vice versa. Since no mechanism fixes the interface at a certain position, 

the interface position is subject to diffusion due to thermal fluctuations. We 

choose the coordinate system such that the origin of the ?/-axis is at the center 

Do of the interfacial profile, which we determine via||51|] 

yo+20 

V] m(y)\ = min. (1) 
yo-20 

Here m = p A — Pb is the order parameter of the demixing transition and the 
relative monomer densities are defined by p A B = 4>a,b/4>i where (ft a and 4>b 
are the volume fractions taken by A and B monomers. The simulation box 
contains 32768 polymers of equal chain length 32. As initial configuration, we 
choose a relaxed configuration of a pure homopolymer interface randomly 
pick 1024 chains with the center of mass at distances of less than ±5 from 
the interface with 5 = 3 or 9, and turn them into copolymers. No effect 
of the choice of 5 on the results has been found. The simulation algorithm 
involves random hopping of randomly chosen monomers by one lattice unit 
with Metropolis probability, but no grandcanonical moves, i.e., the number 
of copolymers remains fixed. After an initial equilibration time of 2.5 • 10 5 
attempted moves per monomer (AMM), the concentration of copolymers in 
the bulk at e = 0.1 is 0.05% ± 0.01% - estimates from grand canonical bulk 



simulations suggest that it should be around 0.04%[52]. We average over 86 
configurations in total, where the data for averaging are taken every 10 4 AMM. 
The area covered by one copolymer can be roughly estimated by ttB? b , where 



■g,b - 



b 2 N/ 12 is the gyration radius of one copolymer block. Hence 1024 



copolymers cover approximately 30% of the total interface area 512 x 512, i.e., 
the system is well in the diluted regime. 

We close this section with a brief comment on the self consistent field calcu- 
lations. In a previous paper, we have compared the properties of homopolymer 
interfaces with the predictions of self consistent field theories for completely 
flexible and semiflexible chains with different chain rigidities ^5j. The quali- 
tative agreement between theory and simulation was over all very good. The 
main quantitative discrepancy was found in the interfacial width - at chain 
length 32, the self consistent field theory underestimates the profile widths by 
at least a factor of 2. This effect could not be explained by capillary waves 
alone, but seemed to be a consequence of the short chain length. At temper- 
atures sufficiently below the critical point, so that critical fluctuations do not 
affect the interface any more, the interfacial width gets already comparable to 
the screening length of the excluded volume, on which length scale chains can 
not be treated as pure random walks. This leads to wrong predictions of the 
size of the interface. However, theory and simulations still agree quantitatively 
for other quantities, e.g. the reduction of the total density at the interface. 

Our present calculations are based on this work. We consider single worm- 
like copolymer chains at an interface of wormlike homopolymers. Chains are 
represented by space curves f(s) with s varying from to 1, and the single 
chain partition function of a copolymer is given by 



where W{(f) is the self consistent field acting on a monomer of type j in a 
homopolymer interface. Each space curve is assigned a statistical weight in 

9 



Z 





the functional integral, V{f(-)} = V{f(-)}Vw {^*( - )} withf53[ 

V w {r(.)}=^H5(u 2 -l)exp{-^L £ ds^fl, (3) 

where a is the fixed monomer length, r\ a dimensionless stiffness parameter, 
u = (Na) the dimensionless tangent vector constrained to unity by the delta 
function, and M the normalization factor. Here we choose rj = 0.5, the value 
which best reproduces bond orientations in pure homopolymer interfaces in the 
bond fluctuation model The self consistent fields Wi for such an interface 
of semiflexible homopolymers have previously been determined numerically in 



Ref. ||45|| , based on the Helfand type free energy functional 

(3F = p j dr{ X $A$B + ttt^(^ + ®b - l) 2 } (4) 

with the total bulk monomer density p, the relative monomer densities $i(r) = 
Pi(r)/p (i = A or B), the Flory Huggins parameter \i an d the compressibility 

K. 

The distribution of copolymer segments is calculated by solving appropriate 
diffusion equations for the end segment distributions Q(f, u; s) and Q + (f, u; s) 
for chain parts of length sN < N, which begin on the A side (Q) or the B 
side (Q + ) of the copolymer^], From those one can calculate the density 
of monomers at position s with orientation u via $(F, s) = Q(f, s)Q + (f, 1 — s) 
or $(r, u) = Q(f, u, s)Q + (r*, —u, 1 — s). The numerical treatment is facilitated 
by expanding the functions Q, Q + in Legendre polynomials and including only 
the three lowest moments ^5| . 
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3 Results 



All results presented here were obtained at e = 0.1 or xN = 17, i.e., the 
system is far from the critical point ((xN) c = 2.4 |48| , [50]) in the strong 
segregation regime. Simulations were also performed at e = 0.05; the results 
are qualitatively the same, but the effects are less marked. 

In the following, lengths are given in units of wssl = bj y/Qx, the interfacial 
width of a homopolymer interface in the mean field strong segregation limit. 
The radius of gyration in these units is R g = 4.2 wssl, and the slab thickness 
D = 37.4 wssl- 

3.1 Density profiles and segment distributions 

Figure 1 shows profiles of A and B monomer densities in systems with and 
without copolymers, profiles of just the A and B blocks of the copolymers, and 
total density profiles. As expected in the diluted regime, the A and B monomer 
density profiles are hardly affected by the presence of the copolymers. The 
interface as a whole is still very similar to a pure homopolymer interface. Due 
to the finite compressibility of the blend, the total density is slightly reduced in 
the interfacial region |44]|. The distributions of copolymer monomers A and B 



agree qualitatively with the experimental results and confirm the simple picture 
presented in the introduction: Monomers of type A are more concentrated in 
the /L-rich phase, Monomers of type B in the B rich-phase. The distribution 
is rather broad, a fair portion of the A monomers sticks into the S-rich phase 
and vice versa. 

Distributions of single chain segments are shown in Figure 2 and 3. For 
homopolymers, one finds a relative enrichment of chain ends p e at the inter- 
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face, whereas the concentration py 2 of middle segments (the 16th and 17th 
monomer) is comparatively low there. The total density profile ph is best re- 
produced by the density profiles pi/ 4 of the segments at one and three fourth 
of the chain (the 8th and 25th monomer). Looking at the higher concentration 
of chain ends at the interface, one is lead to suspect that homopolymers tend 
to form loops with two ends at the interface. However, this is not the case, 
as can be seen from the distribution p e _ e of midpoints between the two ends 
of homopolymers. It is strongly reduced at the center of the interface and 
enhanced at the distance of one gyration radius from there (Figure 2). 

Copolymer segment densities show the inverse trend: The middle segments 
concentrate at the interface, whereas the chain ends stretch out into their 
favorite bulk phase. Hence our results are in qualitative agreement with the 



experimental findings of Russell et al |[25|| . As in the case of homopolymers, 
the total density profiles of A and B monomers are almost identical with the 
distribution of the 8th and 25th monomer, respectively, i.e., the distribution 
of segments in the middle of the A or B block. Note that the peak of the 
distribution of middle segments pi/2 is relatively broad, broader than the radius 
of gyration, hence they are not strongly confined to the interface as assumed 
by the Leibler theory (Figure 3). The results of the self consistent field theory 
are shown in the inset. The fact that the SCFT underestimates the interfacial 
width also leads to quantitative discrepancies in the distribution of copolymer 
segments. However, the qualitative agreement is very good, and in the wings 
of the profile one even reaches quantitative agreement. 
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3.2 Chain and bond orientations 

Next we discuss the orientational properties of the polymers. It is instructive 
to consider separately the orientations of single bonds, of chain segments, and 
of whole chains. Chains with rather weakly oriented single bonds can still be 
strongly oriented as a whole, as found by Miiller et al for homopolymers at a 



homopolymer interface []S3, [45 



The orientation of whole chains involves two different factors. First, chains 
may be oriented without volume changes, i.e., the total gyration radius or 
end-to-end vector remains unaffected by the orientation. A weak orienting 
field is sufficient to bring about orientation of this kind. Second, chains may 
get compressed or stretched in one direction. As Figure 4 illustrates, the lat- 
ter effect dominates close to an interface. The mean squared components of 
the end-to-end vector in directions parallel (x,z) and perpendicular (y) to the 
interface are shown for homopolymers and copolymers. The components par- 
allel to the interface hardly vary throughout the system. Perpendicular to the 
interface, the end-to-end vector of homopolymers is reduced in the interfacial 
region, as is already the case in pure homopolymer systems. Homopolymers 
are thus squeezed perpendicular to the interface, and get effectively oriented 
parallel to the interface. 

Copolymers show the inverse behavior, they stretch in the direction perpen- 
dicular to the interface. As obtained both from self consistent field calculations 
and from the simulations, the effect is very strong for copolymers centered at 
about one to two radii of gyration away from the interface, and much weaker 
for those located in the wings of profile or at the middle of the interface (Fig- 
ure 4). One can picture the latter as consisting of two almost independent 
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homopolymer blocks, which hardly feel the effect of being linked together at 
one end. Indeed, the vector connecting the ends of the single A or B blocks is 
on average hardly oriented (Table 1). Only the blocks centered deep in their 
majority phase (y/wssL ~ 7) stretch perpendicular to the interface, since they 
are pulled towards the interface by the other copolymer end (Figure 5). Note 
that Figure 5 also demonstrates, how minority blocks (A blocks in the B phase, 
B blocks in the A phase) slightly shrink by a factor of approximately 0.8 due 
to the hostile environment. The vector Dab connecting the centers of mass 
of the A and B blocks is on average oriented and strongly stretched in the 
negative y direction (Table 1). 

In sum, single blocks are mostly not oriented at all; the perpendicular orien- 
tation of whole copolymers results from the arrangement of the two constituent 
blocks. Diblock copolymers can be pictured as dumbbells 0] consisting of two 



mildly perturbed homopolymer coils. Similar copolymer shapes have already 
been found by Binder and Fried in simulations of block copolymers in the 
disordered phase, far above the ordering transition |39|. 



Consequently, one would expect that the only region, where the local con- 
formation of a diblock differs significantly from a homopolymer conformation, 
is the region close to the link which connects the two blocks. We study this in 
more detail by looking at the orientation of single bonds b. In order to do so, 
it is useful to define a bond orientation parameter 

= ^ . (5) 

Negative q implies orientation parallel to the interface, positive q perpendicular 
orientation. 

The orientations of homopolymer bonds are shown in Figure 6. Like whole 
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chains, but to a much lesser extent, bonds tend to align themselves parallel 
to the interface. Unlike the end-to-end vector of whole chains, the average 
squared bond length (b 2 ) varies by less than 0.15% throughout the system 
(not shown here), hence the bond length is quasi fixed, the bonds get oriented 
without compression. The effect is strongest in the middle of the chains and 
qualitatively the same, but weaker for end bonds. The bond orientation pa- 
rameter can be compared to the average orientation of the tangent vector u 
predicted by the self consistent field theory g=|(3(u^) — 1). The results show 
the same trend as the simulation data for the middle bonds, but do not predict 
any effect for the end bonds (Figure 6, inset), unlike what is seen in the MC 
simulations. This reflects the difference between the tangent vector at the end 
of a continuous space curve, as assumed by SCFT, and an end bond connecting 
two discrete statistical monomers of a relatively short chain, as it is the case 
in our MC simulation. 

In contrast to homopolymer chains, MC data reveal a rather strong depen- 
dence of the orientation of bonds in copolymers on their position within the 
chain (Figure 7a). The bond which links the two blocks is preferably oriented 
perpendicular to the interface - the stronger, the further away its location 
from the center of the interface. Already the bonds next to the link bond 
show a much smaller effect. The larger the distance along the chain from the 
link bond, the stronger the tendency of parallel alignment at the center of the 
interface becomes, and the further it reaches out into the wings of the profile. 
Hence bonds in the middle and at the end of an A or B block behave very 
much like homopolymer bonds. We note en passant that, compared to the 
other bonds, the link bond is stretched by ~ 4% due to the relative repulsion 
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between the adjacent A and B monomers. 

The results obtained by SCFT show the same trends as the simulation 
data (Figure 7b), and agree qualitatively, but the same remarks as above ap- 
ply. Thus, SCFT does not reproduce the slight tendency of parallel alignment 
which is found in the MC data at the ends of the blocks. The orientational 
parameter of the link bond at the center of the profile is q ~ +0.025 according 
to both MC results and SCFT. However, SCFT underestimates the average 
orientation (averaged over link bonds in the region shown in Figure 7a) by 
~ 30%. At distances of several gyration radii R g from the interface, the bond 
orientation parameter q drops back to zero, indicating that copolymers that 
deep in the bulk are oblivious to the effect of the interface. Note that two 
different length scales are reflected in the profiles of q. One of them, the width 
of the homopolymer interface, fixes the width of the central dip in the profiles, 
the other one, the gyration radius, determines the overall width of the region 
with nonzero q. 

3.3 Self and Mutual Contacts 

Another instructive quantity in the bond fluctuation model is the number of 
monomer "contacts" , i.e., the number of interacting monomer pairs. It is expe- 
dient to distinguish between "self contacts", i.e., contacts between monomers 
belonging to the same chain ("intrachain contacts"), and "interchain contacts" 
of monomers from two different chains. The self contacts provide additional 
information on the conformation of single chains, and the interchain contacts 
on the arrangement of chains relative to each other. 

Figure 8 compares the number of self contacts per monomer N itSe if/pi in 
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homopolymer (i = h) and copolymer (i = c) chains. The main contribution to 
Ni,seifl ' Pi comes from the two direct neighbors of a monomer in the chain, thus 
Ni,seif/Pi is mostly slightly larger than 2. Note however that direct neighbors 
do not always interact with each other, since the maximum bond length ^/Wa 
is larger than the range of the interactions \/6a . The number of self contacts 
in homopolymers is almost constant throughout the system. It is slightly 
enhanced in the wings of the profile and decreases again at the center, where 
A and B meet - suggesting that homopolymer chains tend to loop away from 
the interface, such that there is no room for intrachain contacts at the center, 
but that the number of contacts increases right next to it. In copolymers, 
one finds the opposite behavior - the relative number of intrachain contacts 
is higher at the center of the interface than in homopolymers, but diminishes 
rapidly as one moves away from the interface. This trend is promoted by three 
factors. First, far from the interface copolymers stretch towards the interface, 
and the number of self contacts associated with backfolding goes down. Second, 
the stretching goes along with a higher occupation of the (0, 3, 0) bond, which 
is outside of the interaction region, hence monomers gradually lose contact to 
their direct neighbors. Third, the further away one moves from the interface, 
the greater becomes the contribution of chain end monomers, which only have 
one direct neighbor in the chain. Obviously, it would be highly desirable 
to study separately the contributions to N itSe if/pi from monomers which are 
direct neighbors in a chain, and from monomers which are further apart from 
each other. Unfortunately, such a distinction was not possible due to limited 
memory space. 

Whereas the number of self contacts reflects the conformational properties 
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of a chain, the demixing of A and B chains in the melt is essentially driven 
by the interchain contacts. The number of interchain contacts per monomer 
can be interpreted as an effective coordination number z e ff. In homogeneous 
systems, the identification \ = ^ z eff^ makes contact with the Flory-Huggins 
theory pU[. In inhomogeneous systems, the effective coordination number is 
position dependent. As long as the finite range of interactions is neglected, 
mean field theory simply asserts that it is proportional to the local density of 
monomers: 

Zi,eff{r) = Nl 'p'^ <* ( 6 ) 

A more elaborated mean field approach predicts for a system which is inho- 
mogeneous in one direction y |57[|, 



wW + ^ w with (7) 



where V(r) is the integrable part of the interaction potential {i.e., excluding 
the hard core part) and 7(f) the normalized pair correlation function. Since 
the interaction region extends no further than \y\ < 2ao, an upper bound for 
the factor k 2 is given by k 2 < 4ag. Figure 1 shows that p"{y)/p{y) < COl/ap. 
Hence the expected deviation from the simple proportionality law (|6|) is of 
order < 2%. Furthermore, one would expect that the effective coordination 
number does not depend on whether a monomer belongs to a copolymer or a 
homopolymer. 

In contrast to these considerations, the Monte Carlo data reveal a much 
lower dip of the reduced coordination number z e ff(y)/p(y) (Figure 9) at the 
center of the interface. The discrepancy between mean field assumption and 
simulation data becomes even more manifest when looking at the number 
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of contacts Nab between A and B monomers. From the usual mean field 
assumption, one would expect Nab oc Pa ■ Pb (Figure 10). The simulation 
data however show that Nab/{paPb) is increased by a factor of 1.5, compared 
to the bulk value, at the distance of approximately two radii of gyration from 
the interface, where the chains are stretched (cf. Figure 5). At the center of 
the interface, on the other hand, it is decreased by 25%. 

One can deduce that chains are compactified in the interfacial region. In- 
terestingly, this holds particularly for copolymer chains. We recall that the 
number of intrachain contacts in copolymers is enhanced at the interface, 
whereas the interchain contacts are obviously suppressed. Thus copolymers 
offer relatively little "active surface" for interaction with external monomers in 
the interfacial region. In contrast, the relative number of external contacts in- 
creases further away from the interface. Hence the active surface per monomer 
will presumably increase with the copolymer chain length. This may be an 
additional reason why short copolymers are relatively poor amphiphiles. 

3.4 Summary 

We have studied the ternary system of A and B homopolymers and symmet- 
ric AB diblock copolymers close to an A/B interface in the dilute regime, 
where copolymer coils do almost not overlap with each other, by Monte Carlo 
simulations of the bond fluctuation model. In this regime, the structural prop- 
erties of the interface are not altered with respect to the pure homopolymer 
interface - in particular, we find no significant broadening of the interfacial 
width, as is expected at higher copolymer concentrations. According to our 
simulation results, copolymers resemble oriented dumbells with the A block 
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sticking into the A rich phase, and the B block sticking into the B rich phase. 
The conformations of single blocks are not very different from conformations 
of pure homopolymer coils. Single blocks tend to orient themselves parallel 
to the interface, like homopolymers, whereas copolymers as a whole are ori- 
ented perpendicular to the interface. When looking at the orientations of single 
bonds, we find that the orientational properties of the central bond, which links 
the two blocks, differ distinctly from those of a homopolymer bond, but that 
already the neighboring bonds behave very similar to homopolymer bonds. 
Hence the chain loses the memory of the link very rapidly. Our results are 
reproduced qualitatively, but not quantitatively, by a self consistent field cal- 
culation. Furthermore, we found that copolymers are unusually compact in 
the interfacial region, i.e., monomers have more contacts to monomers of the 
same chain, and fewer contacts to monomers of different chains than monomers 
in the bulk phase. As a result, the copolymer-copolymer interchain contacts 
are particularly suppressed at the interface. This should become important as 
the concentration of copolymers is increased. Future work will be concerned 
with the structure of interfaces with higher content of copolymers and with 
the transition between different regimes, from the dilute to the wet brush to 
the dry brush regime, until the two phase region breaks down and a lamellar 
phase emerges. 
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Table 



Table 1: Orientational properties of copolymer chains. Averages are taken over all 1024 
copolymers and over 86 independent configurations. Compared are x, z components (to- 
gether) with y components. Lengths are given in units of wssl- 
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Figure Captions 



Figure 1: Monomer density profiles as a function of the distance from the 
center of the interface y, in units of wssl = b/y/6x = 1.71 a . Profiles 
are shown for the density of all monomers (p), of A and B monomers 
separately (p^ and p#), of just homopolymer monomers (p/J, and of A and 
B monomers belonging to a copolymer block (pa,c and Pb,c)- Also shown 
for comparison are A and B monomer profiles in a pure homopolymer 
system (p A and p%, taken from pip. 



Figure 2: Homopolymer segment density profiles vs. y/wssi in units of 
wssl = b/^/6x. Profiles are shown for the density of monomers in 
the middle of the chain (pi/2), at the end of the chain (p e ), at one and 
three fourth of the chain (pi/4), and of all homopolymer monomers (p/i). 
Also shown is the distribution of midpoints between the two ends of ho- 
mopolymers, p e _ e . Units of densities are the total bulk concentration of 
monomers p&, or p sb = pb/16 (as indicated). 

Figure 3: Copolymer segment density profiles vs. y/wssL in units of wssl = 
bj \/6x = 1.71ao- Profiles are shown for the density of A and B monomers 
in the middle of the chain (pi/2, squares), at the end of the chain (p e , 
circles), at one and three fourth of the chain (pi/4, diamonds), and of all 
copolymer monomers (p c , broken line). Inset shows results of the self 
consistent field theory for A monomers. Full lines show the predictions 
for segment density profiles pi/ 2 , p e , pi/ 4 , broken line shows total density 
profile p c , and symbols compare with MC results (symbols like above). 
Units of densities are p&/2, or p s b = pb/32 (as indicated). 



Figure 4: Mean square end-to-end vector components (R 2 ) with % = x, y, z in 
units of the average bulk value b 2 N/3, plotted vs. the distance of the cen- 
ter of the end-to-end vector from the interface y in units of wssl = bj \f§X- 
Results are shown for homopolymers and copolymers, and compared to 
the prediction for copolymers of the self consistent field theory. 

Figure 5: Mean square end-to-end vector components (R 2 ) (i = x, z, or y) 
of the copolymer blocks in their minority phase (A block in B phase, 
B block in A phase), and majority phase (A block in A phase, B in B 
phase), in units of the average bulk value b 2 N/6, plotted vs. the distance 
of the center of the end-to-end vector from the interface y in units of 
wssl = b/ y/Ex- 

Figure 6: Orientational order parameter q for end bonds (dashed line) and 
middle bonds (the bonds connecting the 16th and 17th monomer, solid 
line) in homopolymers, vs. y/wssi in units of Wssl = b/y/6x~. Inset 
shows the prediction of the self consistent field theory. 

Figure 7: (a) Orientational order parameter q for end bonds, link bonds (link- 
ing 16th and 17th monomer), bonds next to link bonds (15th to 16th 
monomer and 17th to 18th monomer), and bonds in the middle of a 
block (8th to 9th monomer and 24th to 25th monomer) in copolymers, 
vs. y/wssL in units of wssl = b/y/Ex- 3 ) Prediction for q of the self 
consistent field theory for s = (0,1) (end bonds), s = 0.5 (link bonds), 
s = (0.25,0.75) (block middle bonds). 

Figure 8: Number of intrachain contacts per monomer N iySe i f / Pi(y) (i = h,c) 
vs. y/wssi in units of wssl = b/^/Qx f° r homopolymer and copolymer 



chains. 

Figure 9: Normalized effective coordination number z e ff(y)/p(y) vs. u/wssl 
in units of wssl = b/^/Qx, for homopolymer and copolymer chains. 

Figure 10: Normalized number of AB contacts Nab/ (PA(y)pB(y)), and num- 
ber of AB contacts per monomer Nab/p(u) (inset), vs. y/wssL in units 
of u^ssl = b/ ^[§x. 
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